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iBSTRiCT 


These investigations are aimed at evaluating the 
effect of dielectric substrate on the radiation from 
microstrip antennas. Two approaches are considered* 
ilrstly radiating edges of the microstrip antenna are 
modelled by line and sheet sources of magnetic current* 
Radiation from these magnetic current sources is eva- 
luated by using steepest descent method in Pourier 
transform domain* The second approach is based on 
Weiner-Hopf techniques, A dielectric filled planar 
waveguide with top conductor truncated at 2i=»0 is ana- 
lyzed* Reflection coefficient and radiation pattern are 
calculated* Edge resistance and equivalent extension 
of periphery to account for fringing field is calcu- 
lated therefrom* 
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CHAPTER 1 


INTBDDUCTIOH 

1*1 Mlorostrip Antennas : 

Microstrip antennas are used in vanious comunication 
systems ^ . One basic form of microstrip antennas is shown 
in Pig. 1.1, It consist of a rectangular conducting patch 
in one side of the dielectric substrate backed by a ground 
plane on the other side* Several other configuration as 
circular, triangular, ring type etc are also used^^^. 

To study radiation from microstrip antennas, the rect- 
angular configuration shown in Hg. 1.2 may be considered. 

The cross sectional view of the antenna is shown in lig* 1.5* 
The spacing between the conducting patch and the ground 
plane being only a small fraction of the wavelength, the 
electric field can be assumed to be constant along the hei^t 
of the antenna. The field variation along the length of the 
antenna is shown in PLg, 1.3,* The length is very nearly 
equal to ’^g/2, where Ag is the effective wave length in the 
dielectric medium. The field may vary also along the width 
depending on the position of feed point(and the frequency of 
operation). Radiation takes place from the fringing fields 
which are present at the aperture foimed by the boundaiy of 
the patch conductor and the ground plane# 
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An aperture approach was used by' James and Wilson^^^ 
to approximate the radiation fields of a microstrip antenna. 

In this approach, the total field in a aperture is ' 
considered to be sum of an incoming TEM wave and the reflec- 
ted field at the edge of the radiating patch. If the ref- 
lection coefficient is accurately known, the field distri- 
bution in the aperture can be precisely determined. Even 
in absence of exact aperture field distribution, this appr- 
oach yields a resonably good approximation of the fields, 

[31 

Lo et-aLf- J have mcrdeilrd microstrip antennas as cavities. 
The method based/on the fact that the tangential magnetic 
field along the periphery of the aperture is negligiirly •'Bmall, 
The electric field between the ground plane and the patch 
boundary can be represented as magnetic current sheet. Thus 
the microstrip structure can be considered as a cavity which 
has two parallel conducting plates bounded by a magnetic 
side walls. The magnetic current sheet around the periphery 
and its image on the ground plane is used to compute radiaiion 
fields a 

Oarver and Goffey^"^^ have advanced a similar model 
which also considers a microstrip antenna as a resonant 
cavity* Modal expansion technique is used to find field 
Inside the cavity. The main difference between the two 
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Patch radiator 


Fig. 1.2 Top view of a rectangular microstrip 
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models is in the application of boundary conditions on 
the boundary wall. In this model the boundary wall is 
allowed to have a complex admittance. 

In another mo del^^^ the micixistrip antenna is treated 
as a transmission line resonatar with no tranverse field 
variation. The fields vary along the length which is very 
nearly equal to Ag/2, Bach aperture is considered as mag- 
netic mohopole and the radiation fields are computed from 
the mono pole and its image on the groxmd plane. This 
approach gives reasonable result only for rectangular 
(or square) patch antennas. 

I*resent Investigation ; 

In the methods described above, an approximation is 
made while calculating the radiation fields from the 
computed magnetic current distributions. The approximation 
is to ignore the effect of the dielectric substrate while 
calculating the radiation fields from the , magnetic current 
distribution and its image on the ground plane. The main 
aim of the present study is to investigate the validity 
of this approximation, A single radiating slo% is consi- 
dered for studying the effect of dielectric substrate on 
the radiation from microstrip antennas. Using Elrchoffs- 
Huygen»ii principle the fringing electric field may be exp- 
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ressed as magnetic current distribution. 

In Chapter Two, a radiating aperture is modelled as a 
magnetic current sheet (Hg. 2.2), The current distribution 
being present inside the dielectric substrate, M expression 
for radiation fields for a line current sounce inside the 
dielectric (Hg. 2.1) is derived. Using superposition prin- 
ciple radiation field' from the current sheet is obtained. 

The approach is used also for studying a horizontal current 
sheet placed on the dielectric substrate (Hg, 2,3), 

Chapter Three describes another approach for study of 
microstrip antennas. The edge of a parallel plate wave guide 

I 

where one plate is truncated abruptly (Pig. 3,1) can be con- | 

sidered to be a radiating aperture of a microstrip antenna, f 

TBM wave is assumed to be pragating towards this edge inside 
the wave guide. The calculated reflected field may be used 
to compute edge admittance of the structure. The ¥einer-HQpf 
technique is used for this purpose. The extent of fringing field 
outside the wave guide and outward extension of 
the boundary in the planar resonator model is also computed ! 

from this analysis. 

Chapter Pour summarizes the main results and suggests 
the possible extension of these investigations. 



CHAPTER 2 


RADIiTIOH PROM MAGNETIC CURjSBNT SOURCES OR 
A GROURIED EIEIECTRIO SUBSTRATE 

In this chapter, radiation fields for three types of 
magnetic current sources within (or on) a grounded dielectric 
substrate are derived. The three current sources considered 
are (i) a line current within the dielectric, (ii) a vertical 
current sheet inside the dielectric and (iii) a hoj^zental 
current sheet on the top surface of the substrate. Pourier 
transfoim technique is used for solving the inhomogenous wave 
equations encountered in this analysis. The resultant con- 
tour integration obtained is evaluated by steepest descent 
method for the far field. 

2.1 A Line Current Source Inside the Substrate 

2.1.1 Pormulation of the Problem : 

Consider a dielectric substrate with height b and di- 
electric constant *^'r, with a perfectly conducting ground 
plane on one side as shown in Pig. 2.1 ♦ The substrate 
extends to infinity in both y and z directions. There is a 
line source of magnetics. current at x = a, z = 0, Due to 
latei^l uniformity of the structure, there is no variation 
of field in y direction. Thus = 0 electro- 

magnetic fields have only three non-zero components, 



Fig. 2.1 A line current soure inside the 
substrate . 


X 



Fig. 2. 2 A vertical current sheet inside the 
dielectric substrate • 
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and E’ (H^=H2,j3By=0) , H-fields in regions I and li satisfy 
the following two eq,uations respectively* 

+ iO\ie = 0 (la) 

Hyg =,~^c^eQe^ 6(x-b)6(z) (lb) 

P ^2 >2 

¥here , and 6 denotes Dirac delta 

function. Time variation is taken as . Subscript 1 

and 2 correspond to fields in region I and II respectively. 
The other field components are given by the following 
relation; 



(2a) 

(2b) 


The following Iburier transforms with respect to are 
introduced. 



(3a) 
(3 b) 
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CO 


(x»P) = J 

— oo 

oo 

r 

dz 

(Jo) 

'2^ "1 I 


d^ 

(Ja) 


Taking Fourier transfortr'.s of eqns (la) and il^>) 
two ordinary differential equations are obtained. These 
differential equations are to be solved in the transform 
domain (x,p plane) with the transformed boundary conditions. 
Inverse transform yields the desired field. Taking Fourier 
transform of eqn.(la) and (lb), one gets, 





where 


■ ^o^^l “ ^ X > b 

h^^ y^2 ~ *-5WE:Qey6(x-a) 0^ x^b 

I ■ U l I * n .' 

ilo = f(p^- ^^) 

til = f( Ko^ej. - p^) 

\ = free space propagation constant 


(4a) 

(4b) 

(5a) 

(5b) 

(5c) 


General solution to eqns (4a) and (4b) can .be written 



11 


• +P’e b (6a) 

= B Sin (h^x) + G cos (b^x) a 4 (6b) 

^2 “ -^2 (^ 1 ^) 0 4 x 4 a ( 6 c) 

The set of simutaneous equations given by ( 6 ) are solved 
with the folio vdng boundary conditions. 


Boundary Conditions 

(a) The tangential component of magnetic field is conti- 
noua ‘ at the air-dielectric interface x=b, i,e.' 

i^^(b%z) = 

ITsing ( 3 a) and ( 3 c) one gets 
'■f^(.h*,<x) = VgC'’". 1 ) (7a) 

(b) The tangential component of electric field is conti- 
nous at the air-dielectric interface. 

At x=b, B 2 ^(b'^,z) = S^2^'b~,z) 

Using ( 2 b), (3a) and (3c) this boundaiy condition can 
be written in transformed domain as 





dy 2 


x=»b 


dx 


dx 


(7b) 
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(c) ^ x=a , is continous 


ilo) 


( d) ii x=a. 


Vz 


dx 




Xs;3, 


+ dx 


x=a 




(7d) 


This is obtained by integrating (4b) from x=a"' to a"^. 


(e) Tangential electric field is zero at ground plane, 
at x=G, = 0 


dy 


dx 


= 0 


x=0 


(7e) 


Proper branch of ho should be such that it leads to 
outward propagating or attenuated waves. This re<iuires F in 
(6a) to be ze^ro-. Boundary condition given by (7©) gives A^=0 
in (6c). The set of eq.uations given by (6) can be written 
as 

= De”^^ X b (8a) 


^2 = 2 Sin(h^x) + G Cos (h^x) a ^x ^b (8b) 


^2 ^ A Cos (h^x) 


0 X a 


(8c) 
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Using (8) and the boundary conditions given by ( 7 ), the 
folio vdng set of simultaneous equations (•written in the 
matrix foim) is obtained* 


Sin(h2^a) 

Gos(h^a) 

■^^in( h-j^a) 

0 

■Gos(h^a) 

3 in(h^a) 

Cos(h^a) 

0 

0 

Gos(h^b) 

-Sin(h^b) 

e h 

0 0 

0 

Sin(h2^b) 

Cos(h^b) 

_^-hol> 


J 



H 






0 

0 


D 


0 


Solving eqn. (9) for D, one gets, 


( 9 ) 


-jcoe^e^ Cos(h^a) 
h^ Sin ( 11 ^ 13 ) 


(10) 


Substituting D in (8a) from (lO) and then taking the 
Inverse Fourier Transform, the field in the region x b 
can be obtained. The magnetic field is given by the 
following expression* 

Hyi(x,z) = dp 

- 


(11) 
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where D Is given hy (lO) 

The integration in (Xl) caijBOt' be carried out analy- 
tically, However, the far field can be approximated using 
steepest descent method* Appendix A describe® the steepest 
decent method. The radiation field r ^ Q.) for the magnetic 
line current is given by the expression 


|I"HI' I III ,11 

Hyl ( = "^^o^rW'Tii^ 


e 


J(k^r-V4) 


- Cl2a) 


where fi(») = Cos » Oos (k^ua) 


e^Cos d- Cos (kQub)-ju Sin(k^ub) 


.-fi 


- Sin^ ft) 


(12c) 


©■ and r are cylindrioal coordinates shown in 
Tig. (2.1). 


02.bj 


The radiated power density in free space is proportional 
to the square of the field intensity. Total power (P) radi- 
ated from a unit length of the line source can be written as 
%/2 



(r,a)} 


? ^ 2 


2 


dO- 


(13a) 
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where ri is the intrinsic impedance of free space. 

Radiation resistance is calculated from the knowledge 
of the radiated power for a given strength of the magnetic 
current source. Magnetic current is expressed in volts and 
thus, radiation resistance ’R^' is obtained from the expre- 
ssion 



Pover Radiated in Watts 
Strength of source in Volts 


(13b) 


2,1,2 Numerical Results ; 

Radiation pattern for a line current source is shown 
in Ilg» 2.5. The curves 1,2,3 and 4 are for the following 
cases, (operating frequency 3 GHz). 

1. The source on the ground plane with €^=2.53 and 
b = 0-328 cm. 

2. The source on the ground plane with Cj.=2.53 and 
b = 0,159 cni. 

3. The source on the substrate with = 2.53 and 
b = 0 .159 cm. 

e =1 and the source position from the ground 

'■ T 

plane ^ 0,328 cm. 


4. 




€r = 2.53 
Frequency 3 GHz 


b s 0.159 cm 


b s 0.159 cm 


a = 0,a=b < 0.328 cm 


2.5 Radiation pattern for line current 
source ( IHvI vs 9 ) . 
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!Cal)le 2(a) 

gadlatecL Power from Line Source 

strength of ^.current = X Volt ‘ 
Length of source = 1 m 
frequency = 3 GHz. 


Source on the bottom surface of the substrate (a a O) 


Height of the 
substrate ’b * 
in metre . , 

^r 

Badiated power in 
micro watts 

Radiation 
Besistance in 
ohms /metre 

0^00159 

1.0 

83.26 

12.010 

0.00159 

2.53 

79.63 

12.558 

0.00328 

2.53 

78.10 

12.804 

0.000635 

9.8 

81.566 

12,26 

Source on 

the top surface of the substrate 

Height of the 
substrate *b* 
in metre . , 

^r 

Badiated power in 
micro watts 

Radiation 
Besistance in 
ohms /metre 

0.00159 

1 

82.84 

12.071 

0*00328 

1 

81.49 

12.271 

0.00159 

2.53 

77.80 

12.853 

0.00328 

2.53 

71. 

14.011 

0.000635 

9-8 

80.362 

12.443 
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When is increased the radiation pattern gets modi-* 
fied near end-fire direction. The deviation is more for a 
thicker substrate. By comparing radiation patterns 2 and 
5 we note that radiation patterns for sources on the ground 
plane or on the substrate are <3.uite similar; When the source, 
is on the ground plane the radiation near end-fire direction 
is reduced little more than that of the ' case ■ when 
the source is on top of the substrate* 

The radiated power (for a unit source) and radiation 
resistance are tabulated in Table 2(a). Radiated power 
decreases when dielectric constant or the height of the 

substrate ’b ’ is increased. decrease in radiated power 
indicates an increase in radiation resistance, 

2 A Vertical Current Sheet Xiiside the Substrate ; 

A radiating aperture in a microstrip antenna can be 
modelled as a vertical current sheet inside the substrate 
as shown in Big. 2,2* The radiation characteristics are 
calculated of such a magnetic current sheet are derived in 
this section. 

2.2,1 Bjirmul^ion of the Problem : 

A sheet current can be treated as a superposition 
' •' . of line current sources placed adjacent to each other. 
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Por a vertical current sheet the radiation field, is there- 
fore obtained by integrating the expression (l2a) with res- 
pect to the variable *a’ varying from *0*to*b'. The compo*- 
nent V of radiation field is given by the expression: 


Hy2(r.O): 




(14a) 


Cos & Sin (k ub) 

where f 2 (©)= ™'e^Cos~^ ^osXk^ub) - ju~Sin( k^ub) 

(14b) 


The power radiated can be calculated from the expre- 
ssion (l5a) , Hadiation resistance is calculated from the 
knowledge of radiated power and the strength of total 
magnetic current (in volts) across the sheet. 

2.2*2 Numerical Results : 

Radiation patterns for vertical sheet currents for 
an operating frequency of 3 gHz are shown in Pig,. 2.6 , By 
comparing curves 2 and 3 it may be noted that an increase 
in or in height ’b ' modifies the radiation pattern 
mainly in the end-fire direction* 




1. Er =1 ,b *1/16, U8" 

2. Er =2.53, b = 1/16" 

3. Ef =2.53,b = 1/a" 

4. Er =9,8 ,b sO.025" 








■iSS" 


Fig. 2.6 Radiation pattern for vertical 
current sheets (Hy ys 0 ) , , , 


'f ■' ^ 

' ... . ' ' '6 ' ' 
V y- "f 
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TiBLE 2(b) 

Badiated lEb-ger Prom Vertical Sbeet 


Type of source j Vertical Sheet 
Strength of current ts 1 Volt /m 
Length of source = Im 
Prequency e 3 GHz 


Height of the Radiated Power Radiation Resis« 

substrate b in e in micro watts tance in ohms /m 

metre metre 


•00159 

1.0 

•00159 

2.53 

.00328 

1.0 

.00328 

2.53 

.000635 

1.0 

.000635 

9.8 


.2102 

12.027 

.20 

12.640 

.8895 

12.094 

•8156 

13.190 

.03357 

12.011 

.03273 

12.319 
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todlated power and the radiation resistance are tabula- 
ted in Table 2(b), The current dencity is in -volts/m, and 
the power calculated is for an unit length (one metre) of 
source in y~direction, 

2*3 A Horizontal 6urrent Sheet on the Substrate j 

In a microstrip antenna configuration (shown in Pig. 

^•3 ) there is a horizontal cpmponent of fringing electric 

field. This component can be represented as a horizontal 
sheet of magnetic current on the top surface of the subst- 
rate, Eadiation characteristics of such an aperture source 
are discussed in this section, 

2,3*1 Population of the Problem: 

Horizontal sheet of magnetic current on the substrate 
is treated as a superposition of continonsly placed line 
currents on the substrate. The expression (I2a) which 
gives radiation field for a line source inside the substrate 
is also valid when the line source lies on the substrate. 
The radiation field for two line sources positioned at 
(xsb^ZsO) and x=^b, z- ) is a vector s-um of the fields for 
the two sources individually. It may be noted the two radi- 




Ground piano 


An horizontal current sheet with 
one conducting plane- 


Current sheet 


An horizontal current si 
two conducting planes 







ation fields iia-^-e a phase difference of 2niz Sin ©•/ X^) } 
"being neasured fron thebroad side as shown in jBlg* 

Radiation field for a horizontal current sheet (lig,2,5) 
(top surface without Tpetallization) can be obtained by inte- 
grating expression (l2a) when a=b and z varying from ^0’ to 
'L’j L being the width of the sheet* fhe expression is 

Hyi(r,&)= 

(15a) 

Cos Q. Cos (k ub) 

where f 2 ( 0 -)_ Qqs Oos!k^ub)-ju Sin(kQUb) 



(15b) 




j|Cos(k|^L Sin ©•)-!} 4 Sin (RqI> Sin ©■) 


(15 c) 


2.3.2 Magnetic Current Sheet on a Conducting Surface ; 

A horizontal magnetic current sheet on the substrate 
with top surface metallization [Pig. 2.4] is treated as 
a slot antenna of infinite length in y direction. In this 
configuration, the radiation field fora given current 



distribution is independent of the substrate height and the 
value of The inagnetic current sheet can be expressed 

as a z directed B field distribution on the aperture* 

The maxwell’s equation in transformed domain can be written as 



Solution for ^ in eqn (15 a) oah be written 

~h^x 

Yi = M e 


(16 a) 


(16b) 


with proper branch of h^* ' 

The boundary condition is that the B„ field is zero at 
every point of the top metallized plate except at the aperture* 

At aperture the field is E„ = 1. 

z ■ ■ 

In transformed domain the b*c. is 


--g^Cb, a) = — ^ [e^^ - 1] (16c) 

Prom 16 a and 16 b one can write 

Employing steepest descent the for field may be written as 

i 

^ (9) (16e) 

where f^(&) is given by ( 15 c) . 

2*3*5 Numerical Results : 

Radiation patterns for horizontal magnetic current 
sheet on a dielectric substrate (i) with top surface , 









I able 2(c) 

Eadiated Bower from Horizontal Sheet 
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Type of source: (l) Horizontal Sheet (Top sur- 

face without met alii zation* 

(2) Horizontal sheet with top 
surface metallization 
Strength of current = 1 ?olt /m 
Length of source » im 
frequency = 3 GHz 

( 1 ) 


Height of 
substrate 
metre 

the 

»b» 

^r 

Eadiated Power in 
micro watts 

Eadiation 
Eesistance in 
ohms /metre 

•00159 


1.0 

0.2072 

12.201 

.00159 


2.53 

0.1967 

12.852 

*00328 


1.0 

0.8667 

12,413 

*00328 


2.53 

0,7639 

14.083 

.000635 


1.0 

0.03320 ‘ 

12.145 

.000635 


0,8 

0.32346 

12.465 


( 2 ) 


Width of the 

Eadiated Power 

Eadiation 

substrate *b* 

^r in micro watts 

resistance in 

in metre . . 


ohms /metre 

0*00159 

any value 0*2/0-^ 

12.01 

0.00328 

0.8,94/ 

12,0J 
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netalization (iFig, 2.4 ) and (ii) withoul: metalization (Pig. 

2.3 ) 3-re shown in Pig. 2,7. In the case (i), the radiation 
field in the broadside direction is about half that of with- 
out metalization (case ii) , Por the sheet with no metaliza— 
tion on the top surface (i.e. the case ii) the radiation 
pattern gets modified with increase in e^. 

Radiated power for a current dencity of 1 volt/m for 
both the cases is tabulated in Table 2(c). The length of the 
source is 1 metre. Radiation resistance increases when either 
or height ’b ' is increased, 

2 . 4 Discussion; 

Pour types of magnetic current distributions studied 
in this chapter are: 

(i) A magnetic current line source inside or on the 
dielectric substrate Pig, (2,i). 

(ii) A vertical magnetic current sheet inside the di- 
electric Substrate Pig. (2.2). 

(iii) An horizontal magnetic current sheet on the subs- 
trate without metalization on top surface Pig. (2,3) . 

(iv) in horizontal magnetic current sheet on the subs- 
trate with m.etalization on top surface. Pig. (2.4). 

Radiation patterns get modified by the presence of di- 
electric substrate in all the above cases other than the 
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the type .(iv), when the relative pemitivity of the suhs~ 
trate e^ or height of the substrate 'b^ are increased. Por 
all these cases an increse in or increase in height ’b^ 
reduces the field strength near the end -fire direction. Por 
these three types of source the changes in radiation pattern 
are similar as shown in Pigs. 2.5, 2.8. 

Radiated power also decreases for an increase in or 
in height 'b ' for the sources of the types (i), (ii) and (iii). 
A decrease in radiated power means an increase in radiation 
resistance. Prom the Tables 2(b) and 2(c) it is observed 
that radiation resistance for the horizontal current sheet 
of Pig. 2,3 is marginally greater than that of the vertical 
current sheet as shown in Pig, 2, 2. 


Tables 2(a), 2(b) and 2(c) (for operating frequency of 
3 GHz) show that with £^,=1 the radiation resistance for all 
types of sources l/yvthx.cl^'y^ the type (iv) have magnitudes 
which are in close ' * agreement that obtained f rovi the 
fo xmulaCiT] 

£ ^ .-AS— ohm/m (17) 

1 0 
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where _Xo ^ space wave len^h 

"b = height of the substrate. 

The tables also show that for staall heists with 
®r ^ the radiation resistance calculated has values 
which closely agrees that given by expression (17) • 



CHAPTER 3 


RADIATION EEDM THE. EDGE OE A SEMI-INPINITE 
SHEET ON A' DrELEGTRIC SGBSTRAPB 

Magnetic current soixrce distributions discussed in 
in Chapter Two do not model the micro strip antennas accu- 
rately , Another configuration which can be used to investi- 
gate the radiation characteristics of these antennas is 
that of a semi— infinite conducting plate on the top surface 
of the dielectric substrate. This configuration (Pig* 3,1) 
may be analyzed by using Niener-Hopf technique^ 

3 ♦ 1 AnalrglLo 

3il*l Method of Analysis : 

Consider a parallel plate wave-guide whose top plate 
is terminated abruptly as shown in iig, 3.1. The yz co- 
ordinates are selected parallel to, and x-axis is chosen 
perpendicular to the conducting planes. The ground plate 
extends to whole of the yz plane and the top plate extends 
in y direction from -ooto and in z direction from - oo 
to 0, The wave-guide is filled by a dielectric substrate 
with a dielectric constant of e ♦ let TBM incident field 

■ ■ I, ■ , . X ^ 

be travelling in the z direction as shown in Pig. 3.1. The 
incident field exists in the range 0 4^2^ b and-po ^ 
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and s5ero else where. 

Due to the discontinuity in the top plate, there will 
be scattered fields both inside and outside the dielectric 
substrate. The total field inside the dielectric substrate 
is the sum of the incident cfie Id (exist for ) and 

the scattered field. Outside the substrate as'there is no 
incident field the total field is given only by the scattered 
field. 

One can write = 9^^^ +9 (l) 

where superscript t and i denote respectively total and 
incident field. 9 is the »cattered field. 

Choosing a line variation of e the incident field 

can be written as 



SS8 

with k^> 0, k^> 0| for loss-less case — ». 0*^ and 

02 is conductivity of the substrate. 

Electromagnetic fields are obtained by solving the 
wave eq.uation for soalor potential <p. 
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There is no variation of fields in y directionj thus 
0 and 9 satisfies the follo^ng wave equations. 


+k: ^) 9=0 


0 < X b 


(Vm + k ) 9=0 


X > b 


where Yn 


with ^ 0 and k.^ ^0, for loss lees, case K2 — > ( 
It is assumed that > 1^2. 

various components of the scattered electromagnetic 
field are obtained from the following relations. ¥e have 


'^y ° f = TE 


^ = _=i_ 

wke oke 32 


jj 5. — .. j . — 

z uke ax 

Defining the propagation constant a as a = o + jv, the 
following Fourier Transfoims with respect to z are introduced, 

_ oO 

f(x,a) = y^(x,a) +yjx, a) =~ 9 (x, 2 ) e^“^ dz { 6 a) 
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whe re W ( x, a) = | (p( x, 

y 271: io 


Ztz) e 


daz 


dz 


(6b) 


and Y (x,a) = 

^ 2 % 


'0 


9(x,z)e dz (6c) 


¥> 


oO 

1 r d 


y (x,a)= Y^(x*a) + YJx,a)= ■— \ 4— q)(x, z)e^“®^dz 

i2% J-oo Hx 

(7 a) 


where 


V u,a) =7^ \ ^ 


9(x,z)e^“^ dz (7b) 


•0 


and V ( Xy cc)~ — 9(x,z)e dz ( 7 c) 

J 


oo 


The scattered field 9 attenuates as rapidly as oxp^— k^\zl) 
as. jss,^— ^ CO inside the dielectric. Outside the 

dielectric it is radiated field and behaves asymtoically as 


^1/2 

Where r is the distance from the edge to the obser- 
vation point , 

The scattered field 9 attenuates as rapidly as exp 
(-IC2 1 ^* • Therefore YV(x,a) is analytic 
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in a plane at least for t >--k2 y (x, a) is analytic 

at least for x> Hence y (x,a) is analytic at least 
in the strip < t lc 2 _ 


Equations (3) and (4) are solved by employing Fourier 
transform technique. The boundary conditions are; i) conti- 
nuity of tangential component of total electric and magnetic 
fields across the air-dielectric interface and ii) Zero 
electric field at both the conducting planes. The tangen- 
tial component of incident electric field is zero at x=0 
and x=b. Thus the tangential scattered electric field is 
also continous across the interface. It may be. recalled 
that incident fields exist only in the dielectric region 
0 X b 0 

Taking Fourier transforms of equation (3) and (4), one 

ge t s 



== 0 

0 ^ X ^ b (8a) 

Y^y(x,a) = 0 

x:^ b (8b) 


(9a) 

k^ 

(9b) 
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Following are the 

(8b). 


general solutions of equation (8a) 


and 


V(x,a) = i^^Ca) e 


^(a) 




O^x^b (lOa) 


Y'(x,a) = B-(a) e^^- + B„(a) 


-YX 


x^ b (lOb) 


The Sign of the root of y (branch of y) chosen should 
be such that it leads to an outward attenuating wave. Ihis 
requires that y has positive real part when -k^ < t <k 

The boundary conditions in transformed domain are 


given below 

(a) Tangential component of electric field (total or 
scattered) is zero at the ground plane. 


at X = 0 , < z ^ ^ 

using eqn(s) and (7a) one gets 

yUO, a) 0 ( 11 ^) 

(b) Tangential component of scattered electric field 
is aero at the top conducting plate. 

Ejg w 0| at X =3 b, z <0 
Using eqn(5) and 7(c) 

« Vl (b",a) =0 (lib) 
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(c) langentlal component of scattered electric field 
13 oontinous m the air-dielectric Interface. 


B. 


E. 


^ t at xj='b, 2 > 0 


using eqn(5) and 7(b), 

V4. = V (b*, a) 


(11«) 


(d) Tangential component of total magnetic field is 
continous at air-dielectric interface. 


H, 


(t) 


y 


=fi 

x=b'" y 


, .at x=b, z >0 
x=b » 


Hy(b-,z)+H = H(b+,z) 


Prom eqn(g], 


7 4. j 

Hy(b+, 2 ) _ Hy(b- z) = e 


Taking Pourier transform of above equation. 


y+(»ta) - = — ( 

^ ifzz ->0 


j(k,+a)z 


e 


dz 


lf^(a+k,) 


(lid) 
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(e) Besides the boundary conditions given above, the 
fields have to satisfy edge condition for an unique solution 
of the wave equations. Edge condition requires that the 
electrical and magnetic energy stored in any finite neighbour- 
hood of an edge must be finite. It gives that none of field 
components of (E,H) should grow more rapidly than 
with "c ^0 and r ^ 0 where r is the radial distance 
from the edge Ihe variation of fields near the edge 

can be written as, 

Hy ^ when z — » 0*“ 

Bj, ~ when z — + o'" 

The above two equations give tbs following variations. 


(b,a) when a — » «o 

(lie) 

(b,a) when a — > o© 

Applying boundary condition given by (lla) to (lOa) one gets. 

Also, the proper branch of y gives that in (lOb) 
must be zero. Equations (lOa) and (lOb) are rewritten as 

V(x,a) « A(a) Cosh.(Yi3c) 0 x ^ b (l2a) 

Y(x,a) « B(a) e x ^ b (l2b) 
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Equations (l2a) and (lab) are written in the following 
form at x=b. 

t;.(b-,a;) + Y_Cb- 0.) = A(a) Cosh ( 13 a) 

'h(b''',a) + Y_(b'*',(i) = B(a) e'^'^ ( 3 ^ 31 ,) 

Taking dorlvatiye of (l3a) and (l3b) at z=b, one gets 

( 14 a) 

(14b) 


T jOc) + y__(b ,a) = Sinh (y^b) 
Y+{b\a) +y'(bta) = -TBe "1''’ 


Equations (lib) and ( 14 b) give 

V Sinh .(Y 3 _b) (15a) 

y;(bia) = - YB{a)e -1''= (l5b) 

Erom (llc)^ (I 5 a) and (l 5 b) the following relation 
Is obtained 


-Ey Y B(a)e = Yi^Ca) Sinh. b) (l6) 

Substracting (l 3 a) from (l3b) gives, 

y+(B“,a) + Yjb%a)~yjb“,a) . 

=?B(a)e*"^^ -iCa) Gosh (Yi^) (l7) 


I 5 




. V 70575 



The following relation is defined. 
yjb+ a) YJt~ra) = J-(a) 
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(18) 

Using boundary condition given by (lid) in (17) and (18) 
one gets, 

7==~7“’7”T ^Acc) c B(a) -A(a) Cosh (Yt'*^) (19) 

V 271 ( a+k|^) ^ 

In etjuation (19)# A(a) and B(a) are substitued from 
eqn. {l5a) and (l5b). The substitution gives the following 
relation* 

3 , , ^+(t)'*‘,a) y_i_(b“,a) Cos 

(20) 


Using (llo) the above equation is rewritten as 


V"27i'v a+k^) 


1 , > 1 Co^ ( Y;i ^) 

+ j-(«) = -Y+ {i> ,“) + 7-aI3 TYl'') 


where G(a) 


V ;(»>■,«) ( 

Sihh (y^,^) 







(21a) 


Sinh (Y^b)+ e^Y Cosh (y^^) 

(21b) 



43 


I-et G(a) be factored as (j(a) = a^{a) Ma), 

where a) is analytic in the upper a plane defined by 

T > -k 2 G-^(cx) Is analytic in the lower half a plane 

defined by t ^k 2 . Pactorization of G(a) is described in 
Appendix (B). 


Equation (21a) may be rewritten as 


j(a-k^)Gt(a) 

Hit ( 




Let S(a) « 


j((x— k,) Gr_( a) 

M|wm 

(a + k^) 


a+(a) 


( s^h G-^( a) 
( 21e) 

+ S^(a) (22a) 


where S^(a) is analytic in the upper half a plane 
defined by t ^ — k2 and S^(ct) is analytic in the lower half 
a plan© defined by t < kg. Decomposition of S(a) into sum 
of 04.(0:) and 0 .‘(a) is described in Appendix (c). 


S 4 u'(o:) and 3^(a) are given by the following expressions 


S^( a) ta 


•3 2 k^ 


lf2it (a + k.) 


(22b) 


0 (ce) 


JL 


(a + k,) 


[ (a-k^)(^(a)+2k^ (^(-k^) j (22c) 


I 
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'Jslng (22). (21^,) 

®” following form. 

W_( a) +J.( a) ( a-k^) g .{a)- _ 

(»=+V^r'> <23) 

\ 

is ft; h. 8X1 cl side of* 

half 1 ve eqn Is analytic in the lower 

iialf plane defined hy t < f • 

^ght hand side is analytic 
“ Plane defined by ,x .fc ^ 
unkaowns (j- V *) Jr t>i 2* two 

... ' ' ^ the eqn. Pro:n edge condition (eqn He 

'ind Liouvelia thearem^^^:] ^ 

thearem both sidee of eqn (2Xe) is .ero. 

Thus from right hand side of eqn ( p- 5 ) 

\2j)^ one gets 


y+ (!=-,«)= J.S*’ 2k^ G^k^) 

— (“+v %(“) 


fate 


J e^b 2 (j^(k ) n( j^) 


(24) 


(15a). (15b), (uo) and ( 24 ), 4 {«) and B (n) 
are given by the following relations. 


A( a) 


(3 Cyb 2 (3j|.-(k^) 

/ST'^Sinhy^b 


(25a) 


jb 2 k, <J4,(k,) 

B(a) 88 - ^ c I I 4 , - G.j.(tt) 

/2i Y e*^^ ^ 


(25b) 
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Radiation Pielf^ 


To find scattered field, (25a) and (25b) are substituted 

in (l2a) and (12b), and then inverse Pourier transforms are 

taken. 

In the region x h one gets 


9(x,z) 




B(a) ^ 


da 


(26a) 


as "•23k^b (lj.(k^)-^ 



G^.(a) 



e -(yj' + 


(26b) 

The above integration cannot be evaluated analytically. 
However, the fiir field may be obtained by using steepest 
descent method. Steepest descent method is described in 
appendix ( A) • Employing this method the for field can be 
written as, 

1/2 

Hy s 9( X ,©) a 2bk^ ® ’ ‘ ’ M-icSin ft) 

(26c) 

kr>> 1, X > b 

The cylindrical co-ordinate system used in the above 
equation is shown in Hg (3.1) , 
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Mflectio fi Coeffioient. ; 

Scattered field in region x < b is obtained by taking 
inverse Pourier transform of equation (l2a) . 



(p(x,z)i=: 



00 


f 


-cso 


4(a) CosL (TT3c)e d. 


(27 a) 




S3£. 


-2e bk, 
r ^d 


CJ+(k^) 


1 

2itJ 


J 


Cosh •( x) 
Yi Sinh (y^6) 




4>0 


da 

(28b) 


Using %( a)ssG'( a) /Q;^(a) j eqn (28b) can be rewritten as 

-Oaz 

’ST" 

(28c) 

wfeore u(a) » — (28d) 

Yi Stnh(Y^fe)+er'Y 

Por Z <0, the contour in eqn (28c) may be closed in 

upper half a(a» «4-jz) plane. The integrand has a pole at 

-jk^z 

a*k^# This residue will give a field variation of e 
(Time variation is ). Por loss less case — >■ 

and this variation gives a propagating mode which travels 


9Cx|Z)a-€y k^ Q':f(k^) 



-oo 


Gosh ( Y^x) u( a) e 
( a+k^j ( a-k^j (^f 
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in -z direction. From eqn (28b) the other poles enclosed 
aape sit a=s;|yj^, where 


j " iitoi* 

+ ^{nit/b)2. 


, n = 1,2,3 (29a} 


These poles give attenuated mode of is real. Thus, 
there will be only one propagating mode as long as the follo- 
wing condition is satisfied. 




or 



(29b) 


Assuming, the above condition being satisfied, the field 
inside the parallel plate and for long values of negatix-^ z has 
the variation of e ® • This corresponds to the reflected 
wave. Denoting reflected magnetic field by H (x,2) one can 
write 


H (x,2*) m <» Gt^k^) 


G- (k^) 


(30 a) 


The above expression also can be written in the following 
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Hyj.(x,2)=_e^ ) 2 ^-dV 


(30 b) 

iteflootion ooemolent i. 

^ ^ field B+. Expression (30b) 

XB til© 3*^©f*l0Gt ©d H fiplV? "Pnn^ in • 4 . • 

a Held for unit incident H field. Ibus 


-fc. j.)c;jL.u t i. 

reflection coefficient R ’ r r ' \ 

V it - ~ ^ can be written 

B 


aa 


R’ 


_ 1, , 2 -“JkjZ 

^^r (k^) e i 


■«' “ (30o) 

It is noted timt in ©mrpssinn 4.b. 4.^ 

“ expression C30c) the time variation 

. One can write reflection coefficient In time 
variation by replacing all +j by -j in ( 30 o). Ibua 

«f lection Ooofflolent £' is gi^n by 


Is 


Rs* conjugate of R » (time e^*^*) 


(30d) 


3*1*4 Edge Impedance* 

Oonei daring the trtmcated parallel plate wave guide as 
a transmission line (Pig 5.2), the discontinuity effect due 
to truncation of the top plate can be equivalently expressed 
as a Impedance 2 ^ connected at the plane of the edge (z 5=G), 
Impedance per unit width along x can be written as 

“ ‘=0 Hr ( 31 a) 




f'ig.3.2 Transmission tine equivalent at the edge 
of the truncated parallel plate waveguide 


Fig. 3, 3 Addition of a line section to represent 
edge capacitance . 
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where Z, 


'0 


” t = f 


(31b) 


a = Width Of the wave guide (along y) 

b = heij^t of the wave guide. 

It ie to be noted that reflection coefficient is a funot- 
ion of k^b, and Z^ depends on a/b ratio. 

The impodahoe can be represented as parallel ooBblna 
tion of a resistance and a reactance X^,. Ihe parts Eg 
and Xg can be obtained using the follovdng expressions 



5’b “ % + dBg 







(32) 


Prom ^ or j£j, the edge capacitance (or inductance 
I#g) can bo calculated from the Isnowledge of the operating 
fre(iugncy* fh© following two expressions give C^, n-r 
depending upon the type of reactance. 

** 


( 33 ) 
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Eqiulvalent Extension: 

The reactance (usually capacitive) part of the edge 
impadancG can be represented as a section of open circuited 
transmission line having same characteristic impedance aS . that 
of the parallel plate transmission line. In Pig 3.3 the edge 
capacitance Cj, is shown as a section of the transmission 
lino of length 
3*2 Hmaurlc?-il iiosults i 

TypicifL radiation patterns for radiation from one 
adga of a microstrip antenna modelled as a truncated 
wavoguido (Pig, 3*1) are plotted in Pigs 3*4 and 3.5 fcr some 
values of and e^* Prom these patterns it is observed 
that an incroaee in reduces the field strength near end 
fire direction. In the other opposite direction (©=-90°) 
the field strength increases, Another noteworthy 'feature 
ia the fact that the radiation peak is not in broadside 
(9«0) but around ^40° towards the end-fire direction. 

Variation of Bafloction Coefficients with k^b is shown 
in llg* 3*6, It is seen that reflection coefficient changes 
both in magnitude and phase with i^b, higher the kQb value 
smaller la the rofLeation ooeffieiant in magnitude, 

i'rooi reflection coefficient ’h’ s-ad characteristiG 
impodunco of tho lino, the edge impe’dance is calculated* 

Thy impy dance 2^ is found to be capacitive for the range 
of k^b (,1 to 1) considered. The edge resistance Eg and 
edgy oapacitancG 0^ are plotted in Pigs 3.7 and 3.8 
ruopectively. 
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Edgo capacitaaace 0^, can be expressed as a trans— 
loisoion lins of length, /ikl* The normalized variation of 
oxtonslon i^l/^ with is shown in fig» 5»10# 

In Fig. 3*9 a comparison is made among the valiaes 
of edge or aperture resistance being calculated using 
tho following methods (with operating frequency 3 GHs) * 


(1) by the formula [17] 


E » 
1 


xao A g 

(V) 


24 


(2) by the procedure described in Chapter Two 
(with te 2.55 for vertical sheet current) 

m . ■ ' 

(3) W®iaer*4fepf technique (with e^, = 2.49) 

Tho curves aro numbered accordingly • For lower value 

:\ 

of k b, the results are similar but they do not agree 
for higher value of k^b* Because of the nature of 
analysis procedure, the results based on Weiner-fiopf 
%thod should b© considered most precise* 






Fig. 3, 5 Radiation pattern for the truncated waveguide with 








kS^SSI^-- 1' 
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Wciner Hopf method . 

for Cr = 1.0 0 ;' 


requcncy 3 GHz 
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(1) Microstnp opcncnd formula [' 

(2) By Weiner Hopf method' , "/■ 

(3) Planar microstnp model 




C5HAPTER 4 


a)muDm BsMifiKB 

TMs chpater contains summaiy of the results reported 
In this thesis and some suggestions for further inires- 
tigations* 

4»1 aiMMiRT OP jR||3T]T.TO « 

Ihe Investigations carried out can be devided into 
two parts i) Badiation from a microstrip resonator edge by 
considering equivalent magnetic current sources and 
ii) Study of radiation from a microstrip resonator edge 
based on Veins r-Hopf techniques. 

4*1 #1 jfagnetic Current Sources* 

lo study the radiation characteristics of an edge 
of a microetrip antenna* the following types of equiva*- 
lent magnetic current sources have been considered: 

1) line eoufce located inside the substrate, ii) vertical 
magnetic current sheet inside the substrate, iii) hori- 
zontal current sheet on the top surface of the dielectric 
substrate (Etg*2*3), iv) horizontal sheet on the top 
BjetaUized surface of the substrate (Pig; 2 , 4 ), 

(i) line Current Sources: 

It is observed that for magnetic currents of equal 
Strengths, a line current source located on the ground 
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piano radiates laore power than that for the oase ^en 

the line source is on the top surface of the substrate. 

The percentage difference in powers radiated in these 

two cases is about 2.2 percent for the case' When e «1, 

r ' 

height of substrate 0.328 cm and the frequency is 3 GfHz. 
This difference may be attributed to the fact that when 
the source is on the ground plane » the image with respect 
to the grounded plane produces a field in phase with the 
source field. When the air is replaced by a substrate 
with ej.ta2.53 the difference in the radiated powers is 
about 9 percent. The radiations patterns for these two 
positions of the source (for the thickness of the substrate 

ft 

considered i.s. b .) are very similar. Por £^,=1 the 
3dp points occur at angles 463° from the broadside. Thus, 
the effect of the dielectric substrate is to reduce the 
total power radiated and to increase the directivity in 
broadside direction. 

ii) Sheet Current Sources: 

Comparing a horizontal sheet of magnetic current 
on the top surface of the substrate with a sheet placed 
Vertically inside the substrate,, it is observed that more 
power is radiated in the latter case. The difference is 
about 2.5 percsnt( strengths being equal), Por e^=2.53t 
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this difference increases to about 6.5 percent . 
fhis may be explained by considering the current sheet 
aa combinations of contigously placed line sources. 
for a vertical sheet there are more number of line 
sources which are close to their images than in the 
case of a horizontal current sheet, 

teine3>€opf Appmanh . ■ 

She following noteworthy points emerge from the 
naethod based on Weiner-Hopf technique. 

The radiation pattern for small k^^b is nearly 
flat from - 90 ° to 90° from the broadside (for = 1 
and « 2.49) • Por = 2.49 the peak occurs at an 
angle of around 40 ° from broadside* is k b increases, 
the radiated power in the positive & direction increases 

The edge resistance fig has been defined as the 
resistance corresponding to radiated and surface wave 
power leaving the edge. This resistance decreases with 
an increase in frequency, fhe values of edge resistance 
computed here differ from the calculations based on 
modelling the radiating edge as slot aperture in a 
ground plane [17]. In this slot aperture model the 
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radiation pattern la symmetric with respect to the broad- 
side plane. Employing Woine3>fiopf technique we note 
that the radiation pattern is not symmetric for hi^er 

values of k^bCk^b^ 0.5)* 

... ^ 

ihi.nging field at the edge can he expressed in terms 
of an extension of the physical periphery such that the 
parallel plate capacitance added is equal to the fringing 
field capacitance associated with the edge* This exten- 
sion normalized with the thickness of the substrate has 
boon computed and is found to decrease with increasing 
frequency. The computed results have been compared with 
the following two cases i) the values that obtained by 
using the formula i) = 0*412 (' 5 ' 2g0 ) width 

w 1., and ii) extension associated with the edge 
capacitance of a wide microstrip line (2.^ =5-^ )• The 
extension in a microstrip line also decreases with 
increasing f frequency but they agree closely only -sdaen 
k^b is very small (Lees than 0*1). 

4*2 Suggestions Eor Eurther Investigations ; 

In the equivalent magnetic current approach the 
surface wave excitation has not been considered* If 
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aurfac# wav© power is evaluated, the radiation effi*i- 
onoy can bo computed. Eesistance of a radiating edge can 
be obtained by combining the radiation resistance with 
a resistance which represents the surface wave power, 

Bdge resistance Hg computed by the method based 
on Wein©a>-fiopf technique can be separated into two parts 
one of which is associated with radiated power and the 
other on surface wave power. This can be accomplished 
if radiated or surface power is calculated. 

The edge impedance calculated by Weiner-Hopf method 
may be Incorporated in the analysis of microstrip ante- 
nnas of varioxis shapes, , 



APPENDIX A 
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STEBIBSI DESCENT METHOD 

This appendix describes steepest descent method which 
is a asymptotic integration technique for the type of inte- 
gral given below. . 

■ -h X 

Hy(x,z) = ^ j D(p) e e ° dp, (il) 

— •0 

#,are 

This expression is used in chapter 2, to evaluate magnetic 
field in the air region from the line current source. There 
is a similar expression also in Chapter 3. 


Using the polar coordinates Pig 

X *sar Cos 0-, z c= r Sin 0 - . (A3) 

One can write ( AL) as 
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For analytic convenience it can be assumei that the 
me diutn is slightly lossyf that is has a very small 

positive imaginary part, Ihe branch cuts, sho-wn in Fig 11 
satisfies the condition given by {A 5 ). 

Let, the folio Tfjing transformation be introduced, 

Paek^ Sin =s k^ Sin (5,^ .+ on) (A 6 ) 

and its inverse transformation 

0 o 3in"^(p/kQ) (A7) 

The inverse sin is a multiple valued function and thus 
to ensure one to one mapping from p to 9 plane, the inverse 

transformation is expli city defined as 

6 + h , . 

0 -.3 ln( j ° 

In (18), the logarithmic function is to be intrepreted as 
principal value type with argument of the function lying 
between -nc. to w. 

Transformations given by (IS), represents a mapping of 
the complex p-plane into a connected strip of the complex 
9 -plane. Fig ( 41) and ( 12 ) show the mapping of the two 
sheeted Eieman plane in to 9 — plane. The four q.uadrants of 
the to? and the bottom sheets of the p plane map into the 
regions denoted by T^(top) and p^(bottom) where 1=1,2, 5, 4. 
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The oontooi* of integration path, which is along the real axis 
in p-plano is indicated by in Fig J\l, The following re- 
lation is satisfied for the entire path of 0^* 


h^w f ( k^ Sin Cos 0 

Substituting (-AS) in (M) one gets, 

1 r 

Hy(r,a)a Coe ^BCk^ Sin 0 ) e 


(19) 


( 110 ) 


where g(0) = 3 Cos (0-©') (ill) 

The steepest descent method can be applied to (llO), , 


The idea is to deform the path C^ in (HO) into a new path 

0 , so that the integration can be approximated by integrating 

s * 

Over a Small portion of C^# The criteria on the selection of 
this new path ar© that He g(0) have its maxim urn . value «6 some 
point 9g and decrease rapidly from The point 9^ is known 

as saddle point and is determined by the following equation* 

= 0 ( 112 ) 

0=0a 

The path of steepest descent passes through the saddle 
point and the equation of the path is given by the equation 

(AL3) 



Cos ( 0^ — ©) Gosh i 9 »l- 






1 Complex 0-planc resulting from 
mapping of p plane • 



She main contribution of the integral given by A(lO) 
oomea from a small segment 0^ around jil =? > 8 , On this new 
path indicated by SD? in Pig A2 integration l(lO) can be 
approximated as ' - ' ■ 


Hy(r,ft) 


JL. 

2 % 



Cos 0 D(kQ Sin 0) 


g(0) 

e d 0 

{al4) 


k^r » 


Using tay lor series expansion the above integration can 
>e expanded around and retaining only the leading term, 
ihe asycntotic approximation can be written as 


Hy(r, 0 ') '^( 2 % k.r ^ 


i .1/2 ^ 


[k Cos D(k^Sin &)] 
o o 

' $ / X ' 

im) 



APPENDIX B 


PACTORIZAIION OP G(a) 


This appendix describes the method for the factorization 
of the function oCa) appearing in Chapter 3. fhen factoriza- 
tion is carried out one can write G(a) as 

G(a) « G+(a)Gja) (Bl) 

where G^e(a) and G:;(a) are analytic and free of zeros in the 
upper and lower half a-planes defined by t k 2 and x k 2 
respectively* 


The function G(a) is given below 

( 2 ) 


.(1) 

5 V 

where = 


G(a) * G'^"^^(a) 

Sinh (Yi^) 


(B2) 

(B3) 


and G^ ( a) = 


G^^^(a) is an entire function and it can be expressed 


Y^ Sinh^Yi ) (T^l^) 


(B4) 


in the following infinite product form* 


G^^^(a)*G^^^(G) 


to 


n- 


_a 

iy 


1- )e 


-jab/n7c ^ 


jab /nil 


W 


n 


nsl 


n 
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(l) 

where a « ± simple zeros of Ua) and i® given 

by the relation 


n 


/ 




Prom (B5) one can write, 

(1). . [ m V . . . 


G' 




Il=:l 


ns=:i 


n=l,2,3 

(B6) 

« ^dab/nm 

dy^ 

(B7) 

-jab /nit 
) e 

(B8) 


with the asymptotic behaviour of G+(a) given as [10] 

( l) /2 

G- (a) ^ ia exp(jba/7t) exp[(-jba/7i) l-l>(-3ba/m) ] 

(B9) 

as a - — ^ 

To factorize the function G^^^(o:), the formula derived 
by Bates and Mittra [13] , is It is given as, 

let G^^^(a) be an analytic function of a, and it satis- 
fies the following conditions in the strip (t| < . 

(a) G^^^(a) is analytic in the strip 

■<( 2 ) 4 cs non ^ero and even, nC2) 


(b) G^ ( a) is non ^s^ro and even, that is G ( -a) — 

.( 2 ) 

( 2) / _^ I “ Ur cf , where a' and h 


( c) G^ ( a) Ba^ e i “ i as 
real ■ 


a 


72 


Then for a within the strip, one can write 


(BIO) 


where Gr^ '^(a) and G ^ (o.) are non zero and analytic in the 


upper (r and lower {% ^ x a-plane, respectively, 

The expressions for G.^^^(a) and a) are 


Si ^\o()=a-^®^(o!)=/9^^ho) (i+S^°^ F 


(1+ f T7 (1+ -f-) X 

m=l m 


llj-(lt j-) exp [-^ + in ( , -Ki(a)+ B^( a)3 

^ (Bll) 


In the above .expressions, the following terms should be 
interpreted as given below 


q(a) = l5:( w)ln [1 + 




where k(w ) = ^ “ [B(w ) + B(~w ) ] 
where B( w ) = -,,ln (G w^) 


(B12) 


(B13) 


(Bi4) 


= Lim wB ( w) 


(B15) 


Zm w Simple zeros of G'^^Ca), where Im 

m = 1, 2| 3 * M* 



13 


Bn 


Simples poles of where Im Pn.„ 




n ss 1, 2, 3 . #« H. 

Rn(a) £ siftn of the residue contributions from the poles 
of B(w ) on the positive real axis. 


/ 


2 2 

Y as Y 0! - k with proper branch 

k « + jk 2 t k^> 0 and k^ 

The asymptotic behaviour of (T^^^(a) is given as 

(B16) 

The branch singularities of are at a = + k , 

( 2 ) 

Prom condition (c) the asymptotic behaviour of o'* ^(a) 
can be written as 

(BIT) 


. ( 2 )/ \ /2 ^ r ibg / 2a n-\ 

Qt+^ '{a) a ' exp [- In { )j 

,( 2 ) 




la )- 


" , Itl C + 
( 2 ) 


Hence V = 0 and h=b. The function '(cc) does not have 
any zeros, but has poles at a = + Pn» ^ ~ N. 

The poles Pa are determined from the equation 


SirbCy^^b) + e^Y Co£^MYib)| = 0 

a = Bn 


(B18) 


Using (BI 4 ) the function B( w) is computed and obtained 


as 


B( w ) 


1 . X + .jy 

u 


(BI 9 ) 
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■■ O' O r\ 

where X «w [h ( w /w») Sinw’b Cos w'b + 1 + Gos^ w'b 

(e^^~1)] {B20) 

y sa Ey Sin. w»b Cos w’b Cos w'b ( ^^/ w* -l)+e^ 

(B21) 

u = (cy wCoB w*b)^ + ( vf' Sin w*b)^ (B22) 

with w’« ( (B25) 

BCw ) does not have poles on the positive real axis) 
and thus Bj^Ca) =0, 

Using (B15) and (B18) one gets 

Vq * 1 (B24) 

Using (B13) and (B18) the function k(w)' is obtained as 

2 

V Sin w»b Cos w'b (2 + e w^b 

K(w) = ^ - ... r— — r (b25) 

^ n (e^w Cos w'b)^+(w* Sin w'b) 

Substituting all the above results in (Bll) the fmction 

/ p\ 

C+' '(a) is computed to obtain 

( a) :=G J ( a) ^ (O) ( 14^ ) ]^(l+— )“^ 

^ n=l ?n 

e xp[- In ( ~^) + q( a) ] (B26) 

where Bh is obtained from (B18)^ and q(a) is to be computed 
using (B25) and (B12). 
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Iho asymptotic behaviour of is obtained by 

pflin^i: (Bl6), The result is, as a- — and x y^xtf: 

^+^^(a) ^ exp [- In ( ^) (B27) 

The function G+(a) can be written as 

C- ■ (J(a) = 0?^“^ (a) (B28) , 

where G+^^^(c£) and are given by (B?) and (B26) 

respectively* The term e is necessary to have an 

unicjue factorization of G(a), The function X(a) is analytic 
in whole of the a-plane. To deter-nine (a) one has to 
invoke edge condition* Edge condition imposes that G^( a) 
and Qi;( a) should have algebric behaviour for large a in 
the respective plane of regularity. 

Using (B9) and (B27) asymtotic behaviour of G+(a) 
given by (B28), can be written as 

, , ' " , i ' 

G-^fia) AjoT^^^ exp [l-G+ln('|2) + j 1) + j 2 j4.^(oc)^ 

(B29) 

ETOm (B29) the function X(a) can be choosen to have 
X(a) = -2^ [ 1 - G + In (||-) + j I 3 (B30) 

where C = 57721 is Euler’s constant. 
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Appendix G 
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r 

DiSOOMBDSI'IION OP PIMCTION S(a) 

I'hio appendix describes the aiethod of decomposition 
nf function S(a) into the forai S(a)=S ,(a) +S (a) where 

'T 

are analytic in the uppef arxd lower half 
oc—plano defined by t “ ^2 ^ respectively# 

Thu function S(a) (Chapter 3) is given by the ex- 
pression, 


S(a) ss 


j(a - k^) G^(a) 
fz% (a + k^) 


do) 


Tht function G (a) is anaiytio in the lower half 

•til 

a-plan© and so only discrete pole at the lower by half 
a-plane is at cxa-k^^* Poliowing the standard method [10], 
[11], one can write, 


S^(a) « 


-j2k^ Gj-k^) 

/2i(a + k^) 


^ 2 k^ G^ ( k^) 

)[Z% (a k^) 

and S^(®) ** S(a) *• S_^(a) 

„ ^ -[{a-ta) Gjo) 9J-V? 

'fz%i<x+ k^) 
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